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Abstract 

The quark-antiquark bound states are discussed using the relativistic specta- 
tor (Gross) equations. A relativistic covariant framework for analyzing con- 
fined bound states is developed. The relativistic linear potential developed 
in an earlier work is proven to give vanishing mesons q + q decay ampli- 
tudes, as required by confinement. The regularization of the singularities in 
the linear potential that are associated with nonzero energy transfers (i.e. 
q 2 = 0, q^ ^ 0) is improved. Quark mass functions that build chiral sym- 
metry into the theory and explain the connection between the current quark 
and constituent quark masses are introduced. The formalism is applied to the 
description of pions and kaons with reasonable results. 
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I. INTRODUCTION 



Description of simple hadrons in terms of quark-gluon degrees of freedom has long been 
an active area in physics. With the advent of Jefferson Laboratory, which operates at inter- 
mediate energies and therefore probes the structure of hadrons, there are new opportunities 
to test simple theoretical descriptions of quark interactions. The first natural step in this 
direction is a thorough understanding of how to treat the relativistic quark-anti quark bound 
state problem. In this context, NJL inspired models have gained popularity in recent years 
HIIJ]. The common goal of these works is to bridge the gap between nonrelativistic quark 
models and more rigorous approaches, such as lattice gauge theory or Feynman-Schwinger 
calculations. While the Euclidean metric based calculations are increasingly popular, their 
applicability, because of the extrapolations involved, is only limited to light bound states 
such as the pion and kaon. Therefore, it is important to develop Minkowski metric based 
models which can be used over a wider scale of energies. One such work using the specta- 
tor formalism was developed in Ref. In those works a relativistic generalization of the 
linear potential was developed and the pion was shown to be massless in the chiral limit. 
However, the calculations involved some approximations and related conceptual problems. 
In this work we improve and simplify the model presented in those works and address in 
detail some of the conceptual issues related to confinement. 

If a quark-antiquark pair (referred to collectively as "quarks") is confined to a meson 
bound state with mass /i, then the bound state can not decay into two free quarks, even if 
the sum of the quark masses is less than the bound state mass. This trivial statement can 
be realized by two possible mechanisms: either (a) the quark propagators are free of timelike 
mass poles, |2| or (b) the vertex function of the bound state vanishes when both quarks are 
on-shell. In this work we prove that the Gross equation supports the second mechanism 
of confinement. The first mechanism, which is commonly used in Euclidean metric based 
calculations, is a stronger constraint since it forbids any free quark states. On the other 
hand, the Gross equation allows one of the two quarks in a meson to be on-shell, but 
insures that the matrix element which couples the bound state to two free quarks vanishes. 
The spectator formalism facilitates the use of the Minkowski metric, and the confinement 
mechanism of this approach has a closer resemblance to nonrelativistic models. 

The organization of the paper is as follows: In Sec. |TJ we review the formalism for 
nonrelativistic confinement in momentum space. In Sec. Ill we outline the general philosophy 
of the spectator approach to the treatment of confined systems, examine the implications 
of confinement for the scattering amplitude, and prove that the relativistic linear potential 
used in earlier works automatically insures that fi —>■ q + q vanishes at the momentum where 
decay of the state into two physical quarks would otherwise be kinematically possible. The 
treatment is first presented for scalar particles, and then generalized to fermions. In Sec. [IV| 
we construct quark mass functions that have the correct chiral limit and preserve asymptotic 
freedom. Our numerical results for pseudoscalar bound states are presented in Sec. [V], and 
some conclusions are given in Sec. [VJ. 
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FIG. 1. The linear potential in coordinate space for e = 0.1 and o = 0.2. The solid line is 
Vs(r), the dashed line is Vz,(r), the dotted line is Vk( r )i and the dot-dashed line is V(r). For 
"small" r < 1/e (the region inside the small box) Vj,(r) and Vs(r) are both approximately equal 
to or. 



II. NONRELATIVISTIC CONFINEMENT IN MOMENTUM SPACE 



We start by reviewing the discussion of confinement within the context of the nonrela- 
tivistic Schrodinger equation given in Ref. [I]. We will denote potentials in coordinate space 
by V and in momentum space by V. The nonrelativistic linear potential is 

V(r) = or. (2.1) 

This potential can be constructed from familiar Yukawa-like potentials in two different ways: 

( Vg(r) = ore~ er (a) 

V ^ = H V L {r) es -°- (e- - 1) = V A (r) + - • (b) (2 ' 2) 

k e e 

These various potentials are shown in Fig. 1 for the illustrative case of e = 0.1 and o = 0.2. 

Note that the two potentials Vg(r) and V^(r) both approximate the linear potential V(r) 
when r << 1/e, but that these two approximate potentials behave very differently at large r. 
The potential Vg(r) — > at large r, so that, strictly speaking, it does not confine particles 
at all. This potential always permits scattering, although when e is small the scattering 
is strongly resonant, and the wave function is significant at small r only for energies near 
one of the allowed resonances. The width of these resonance states becomes narrower, and 
their wave function approaches that of a bound state, as e —>■ 0. In contrast, the potential 
Vl(t) — > 1/e as r —>■ oo and therefore binds particles with energies E < 1/e. Ase^O this 
potential does not permit scattering; it has a spectrum of bound states only. 

Yet for sufficiently small e, it should be possible to move freely from one of these potentials 
to the other, and the results obtained with either form should be equivalent. We will return 
to this later in this section. Now we follow Ref. JTJ and work with V L given in Eq. (|2.2|b). 
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The momentum space form of this potential can be written 



Vl(q) = lim 



VA(q)-5 3 (g) / d^V A {c() 



where 



VA(q) 



87TO" 



2 _L ^2 ' 



Note that the second term (the "subtraction term" ) insures that 

Jd 3 qV L (q)=0, 



(2.3) 



(2.4) 



(2.5) 



which is the momentum space form of the statement that V(r = 0) = 0. The Fourier 
transform of Va is, for finite e, 



Va(t) 



d 3 q 



-iq-r 



Va{<1) 



-a- 



~ lima r 



(2.6) 
(2.7) 



and the subtraction term cancels the singular 1/e term insuring that the linear part of the 
potential has the correct behavior in the limit as e — > and that it vanishes at the origin 
(r = 0). Now, adding a constant potential Vc 



V c (r) 



-c 

-{27i) 3 5 3 {q)C . 



(2.8) 



to the linear potential ( |2.3| ), and inserting the total potential into the momentum space 
Schrodinger equation gives 



2m i 



E 



— Vk(p-k) [*(k,p )-*(p,Po)]+C*(p,p ), (2.9) 



(2tt 



where mjj is the reduced mass, E is the binding energy, and po is an eigenvalue given by 



p = 2m jR £' . 



(2.10) 



The constant potential is used to adjust the energy scale. 

While Eq. ( [2.9|) was derived for the linear potential with the specific choice of Va given 
in Eq. ( |2.4j ), it is instructive to consider it in its most general form where Va is an arbitrary 
function. From this point of view, the role of the second term in square brackets in Eq. ( |2.9|) 
(which arises from the subtraction term), is to insure that the coordinate space potential 
Va{t) is redefined so that it is zero at the origin; ie. Eq. ( |2.9| ) is a standard Schrodinger 
equation for the potential 



V L {r) = V A {r) - V A (0) . 



(2.11) 
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Looking at it this way, we see that any potential Va{t) for which Va{t ) — ^a(O) = oo, for 
some r OJ gives a confined system when used with Eq. i \2.Q ). For example, even the choice of 
a pure Coulomb-type interaction for Va, 



V A (r) 



(2.12) 



would give confinement. The subtraction term forces the interaction to vanish at the origin, 
which requires an infinite shift in the energy (just as in the case of the linear interaction) 
forcing the interaction to go to infinity at large distances. The role of the subtraction is an 
essential part of introducing confinement. This trivial point is worth emphasizing because 
when we arrive at the relativistic equation, the subtraction term will prove to be just as 
crucial as it was in the nonrelativistic Schrodinger equation. 

We know that Eq. ( |2.9|) confines the quarks because it was derived from a coordinate 
space equation which confines, but it is instructive to see in a simple, direct way how 
confinement can be demonstrated directly from the momentum space equation. To see this, 
first consider the case when C = 0, let e be small but nonzero, and rewrite the Schrodinger 
equation 



2m 



-E-V A (0) 



*a(p,Po) 



d 3 k 



VA(p-k) ^a(Kpo) 



where, for the linear potential introduced above 

V A (0) = -- 



(2.13) 



(2.14) 



[For simplicity, we will sometimes refer below to Eq. ( |2.13| ) as the bound state form of the 
equation.] In coordinate space, the potential V^(r) approaches zero at large r, as illustrated 
in Fig. [1]. Hence scattering will take place only if the l.h.s. of this equation has a non-trivial 
solution, which requires 



p 2 >p 2 + 2m R V A (0) =p\ 



Note that this implies that 



E = ^>° 



2m 



oo 



R 



(2.15) 



(2.16) 



as e — > 0, showing that no scattering can take place for finite energies. At energies below 
1/e, only bound states can occur. This is the demonstration we seek. 

Even though Eq. ( |2.13| ) shows that there is no scattering when e — ^ 0, it is still instructive 
to write a scattering equation for finite e. To this end it is convenient to replace Vi(r) by its 
counterpart, Vs(r) defined in Eq. (|2.2| a). This potential has no subtraction, so its momentum 
space Schrodinger equation is simply 



2m R 



E 



^s(P,Po) 



d 3 k 
(2^ 



V s (p-k)^ s (k,p ) 



(2-17) 
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FIG. 2. The potentials Vs(r) (upper solid line), Vl(t) (dashed line), or (dot-dashed line), and 
Vo(r) (lower solid line) in coordinate space for e = 0.1 and a = 0.2. 



This will be referred to as the scattering form of the equation. 

As stated above, we will assume that the two equations ( |2.1ci| ) and ( |2.17| ) give equivalent 
results when e is very small. Their equivalence is clear on physical grounds, since there is 
very little difference, on a sub-atomic scale, between a barrier which is a mile thick and 
one which is infinitely thick. To emphasize this point, Fig. ^ compares the short distance 
behavior of the potentials Vs, Vl, and 

V (r) = V L (r)-V s (r). (2.18) 

As e — > for a fixed range of r, Vq — > and Vs — > Vl- However, a careful mathematical 
treatment of how these two equations approach the limit as e — > presents some subtle 
issues fHH] which we defer to a subsequent paper. Our arguments in the remainder of this 
section are based on simple physical considerations. 

In connection with the scattering form ( 2.17]) we introduce a scattering state wave func- 
tion defined by 

MP,Po) = ^frj^P-p') _ 2m*M a (p,p') (JUg) 

P ~Po 

where Ms is the half off-shell scattering amplitude, and p' 2 = p 2 ,. The wave function 
( P-19| ) has the form of the usual scattering wave function, with the 5 function describing the 
asymptotic plane- wave part. We have chosen to multiply this plane- wave part by a (small) 
parameter rj. This parameter can be removed by dividing the wave function and the half off- 
shell scattering amplitude by rj, so it is, strictly speaking, an arbitrary scale factor. However, 
if we wish to compare the scattering solutions to ([2. 171) with the bound state solutions to 



fl2.13| ), it is necessary to choose rj so that the wave functions are comparable at small r, 
as illustrated in Fig. ^|, and this will require that rj be very small. Such a comparison is 
only possible at certain energies (close to the bound state energies) where the scattering 
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FIG. 3. Comparison of possible wave functions ^>a{ t ) (falling dotted line) and ^s( r ) (heavy 
solid line). [For reference, the potentials Vs(r) (thin solid line) and Vl{v) (rising dotted line) are 
also shown.] The normalization is chosen so that a *-* ^5, making the plane wave tail of ^fg 
(shown in the box) small. In this example rj ~ 0.05. 



solutions are resonant and therefore much larger at small r than at large r. In general, at 
other energies, Eq. Q2.17| ) will have nonresonant solutions that can not be large at small r. 
Only the resonant solutions of the scattering form ( [2.17|) will converge to the bound state 
solutions to Eq. (|2.13|) , and for these r\ is very small. The non- resonant solutions to the 
scattering form are confined to the large r region, and move off to infinity as e — > 0. This 
complicated limiting process will be summarized by the equation 



^s(p,Po) <-> ^a(p,Po) 



[2.20) 



where the <-> symbol means that the spectrum of resonance scattering states obtained from 
( [2.17] ) converge to the bound states obtained from ( [2.13j ), and the nonresonant solutions to 
(p.17 ) can be ignored because they contribute only at infinite energy. 

With this insight, we substitute the scattering wave function ( |2.19| ) into the Schrodinger 
equation ( |2.13| ), giving 



M s {p, p)=V Vs(P - P ) - 2m R j j^ v s{P ~ k) 



Po 



(2.21) 



Alternatively, we may work directly with the bound state form ( |2.13| ) of the equation. In 
this case we make the replacement 



^a(p,Pq) 



2m R M A (p,p') 



P J 



Po 



and, substituting this into Eq. (|2.13|) , obtain the following equation 



J2^ Va( - P ~ k) 



(2tt)' 



M A (k,p') M A (p,p') 



Pi 



P' 



Po 



[2.22) 



(2.2?,) 
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Following the argument developed above, in the limit e — > (and 77 — > 0) the two amplitudes 
Ms and M A should be equivalent. In the notation of Eq. (|2.20|) 



Ms(p,p')~M A (p,p'). 



^2.24) 



We will find it convenient to use Ms when e is very small but nonzero, and to use Ma when 
we want exact confinement (e = 0). Only Eq. ( |2.23j ) has a well defined mathematical limit 
when e — > 0. In our subsequent development we will assume that either Ms or Ma may be 
used with equivalent results. 

When e = 0, the inhomogeneous term vanishes and there exist bound states only. We 
introduce the vertex function 7 defined by 



^(P,Po) 



P 2 -Po 



[2.25) 



The Schrodinger equation for the vertex function, restoring the constant interaction term, 
is 



l(p,Po) = -2m R J ——V A (p - k) 



(2*y 



7(k,p ) _ 7(P.Po) 
,k 2 -Po P 2 -Po 



+ 



2m R C 7(p,po) 



Po 



Next, look at this equation when p 2 — > p\. To this end first write 

7(P> Po) = l(Po, Po) + (P 2 - Po) ^(P> Po) , 
and then substitute this into Eq. ( |2.26|) [with C = for the moment], giving 



^2.26) 



^2.27) 



7(P,Po) 



J^y V A(P ~ k ) 

d 3 k 



k 2 - pg p 2 - pi 



——V A (p-k) [K(k,p Q )-K(p,p Q )} . 



(2.28) 



All terms on the r.h.s. of this equation should be regular as p 2 — > p\. Because of the 
subtraction, the term involving 1Z is finite, and, because of our choice of po, only one of the 
two remaining terms is zero if e is finite 



lim 



lim 

p 2 -*pg 



d 3 k Va(p — k) a f e±2ip 



(2tt) 3 k 2 -^ 
d 3 A; Vk(p-k) 



>o 



finite 



(2vr) 3 p 2 



Po 



00 . 



Hence the subtraction term will be singular unless 

7(Po,Po) = 0. 



^2.29) 



^2.30) 



This condition also insures that the constant term is not singular. We will discuss the 
physical interpretation of this result in the next section. 
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FIG. 4. One of the two-channel Gross equations for the bound state vertex function T. In this 
figure the x means that the particle is on the mass-shell. 

III. CONFINEMENT IN THE SPECTATOR FORMALISM 

A. Introduction 

At this point it is very tempting to generalize the nonrelativistic linear potential Eq. ( |2.3| ) 
by simply replacing the three vector q by a four vector q 

V(q) ± \im[V A (q) - 5 4 (q) f d 4 q'V A (q')} + (2n) 3 5 4 (q)C. (3.1) 

This, seemingly obvious, generalization will not reduce to the correct nonrelativistic limit 
because of the unconstrained behavior of the / dq' V^(g') integral. Lacking a four dimensional 
expression for the linear interaction that reduces to the correct nonrelativistic limit, we 
rephrase our question: Can one find a covariant equation that reduces to the Schrodinger 
Eq. ( p. 26 ) with a linear interaction? The confining relativistic bound state equation should 



be a relativistic generalization of Eq. (12.26 ) . 



A covariant equation with the correct nonrelativistic limit is the Gross equation [§|H. If 
the two quarks have unequal masses m x > m 2 , the one channel equation may be used. It 
has the feature that the four dimensional loop integrals are constrained so that the heavier 
constituent (with mass mi in this example) is restricted to its positive-energy mass shell 
(provided Mb > 0; see Ref. |7|]). However, if the particles have equal mass (mi = mi — m) 
and the mass Mb of the bound state is comparable to m, a symmetrized two channel equation 
should be used. This is illustrated in Fig f|. In this case an average of the contributions 
in which either particle 1 (channel 1) or particle 2 (channel 2) are on their positive-energy 
mass-shell are included, and this leads to a set of equations in which the two channels are 
coupled. The symmetrized two-channel equation has been used previously to describe of 
low energy NN scattering ||. Finally, if the masses are identical and the bound state 
mass is very small (ie. Mb << m), as in the chiral limit, then a four channel equation 
is needed. The four channel equation is a symmetrized version of the unsymmetrized two 
channel equation used in Ref. 0. One of the purposes of this paper is to improve on this 
previous work. 



B. One channel scattering equations for scalar quarks 

We will begin with the one channel equation. The momentum and mass of the quark are 
Pi and mi, the momentum and mass of the antiquark are P2 and m.2, the total momentum 
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is P, and the relative momentum is p, where 



P 

P 



Pi 
1 



+ p 2 = {M B ,0} 
(Pi - P2) ■ 



(3.2) 



The quark will be on mass-shell, and the symbol pf will be used to denote the particle on its 

Ei(p) ■■ 



positive energy mass-shell, (ie. pf = m\ and (pf)o = Ei(p) = yrnf + p 2 ). The scattering 
amplitude M.(pf,p 2 ,p'i~,p 2 ) is denoted Mu(p, p', P), or in the one channel case where there 
can be no confusion, simply by M.(p, p',P). Then, introducing a relativistic generalization 
of the potential Vs, the one channel equation for the scattering of scalar "quarks" (mi > 1712) 
can be written 



M s (p,p',P)=vV s (p,p',P) 



2mim 2 r d 3 k V s (p, k, P) M S (K p', P) 



(2tt) 3 
2m 2 CM s {p,p',p) 



ml 



[P-kt 



ml 



(P-p+ 



(3.3) 



This equation is the relativistic generalization of Eq. ( 2.21 



Alternatively, the bound state form of the scattering equation is 



M A {p,p',P) 



2m\m 2 f d 3 k 



(2n) 



E±(k) 



V A (P,KP) 



M A (k,p',P) M A (p,p',P) 



ml -(P-kt) 2 ml-(P-pt) 2 



2m 2 CM A {p,p',P) 



m 2 2 



P 



Pi 



(3.4) 



This is the analogue of Eq. ( |2.23| ) and has a smooth limit as e —> 0. The kernels Vs and V A 
will be specified later (see Eqs. ( [3.16|) and ( |3.17|) below). Equations (|3.3|) and (|3.4j) will be 
our starting points for this section. 



C. One channel bound state equation for scalar quarks 

In the vicinity of a bound state of mass Mb, or a very narrow resonance with mass and 
width M B = M R + iMj, the scattering amplitude has the form 

ka t ' m T x (p,M B )T x (p',M B ) 

Mx{p,P,P) = M2 _ p2 +Kx(p,M B ), (3.5) 

where X = A or S, depending which of the two forms ( |3.3|) or (|3.4j ) we are using. If e is 
finite and we are using Eq. (|3.3j), the width Mj 7^ 0. If we use Eq. ( |3.3| ) the width is zero 
for all states with mass below some critical mass M e — > 00 as e — > 0. 

Substituting the form ( |3.5| ) into either Eq. ( ^.3| ) or Eq. ( |3.4j ), and equating residues at 
the pole (real or complex) gives the bound state equations for the vertex functions Tx- 
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r s (p, M B ) = -2m x m 2 J 



d 3 kV s (p,k,M B ) T s (k,M B ) 

(27r) 3 E 1 (k) ml-{M B -ktf 
2m 2 CT s ( Pj M B ) 



T A (p,M B ) = -2m im2 J 



ml - {M B - pf) 2 ' 
d 3 k V A (p,k,M B ] 



(3.6) 



(2tt) 3 £i (k) 
2m 2 CT A (p,M B ) 
m%-(M B -pt¥ 



T A (k,M B ) 



r A (p,M B ) 



m\ - (M B - kf) 2 m\ - (M B - p\ 



(3.7) 



where we use a mixed notation with M B denoting both the mass and the four vector {Mb, 0}, 
the difference being clear from the context. 

As with the scattering amplitudes, the two vertex functions are equivalent in the limit 
e -> 



T s (p,M B ) ^r A (p,M B ), 
but the vertex function T A is more convenient to calculate in the limit e — > 0. 



(3.8) 



D. Normalization condition 



The bound state equation and the normalization condition for the bound state wave 
function can be derived from a nonlinear form of Eq. ( |3.3| ) . In this paper the derivative 
dVs/dP^ = in the rest frame, so the result is 

op, d r d*k [ r s (k, Mb)Ts(KMb) \ 

~ dPj (2vr) 3 E 1 (A;) \ m 2 - (P - kf) 2 J ' 1 J 



In view of the relation ( p.8|) this relation can also be written 

2pM d f d 3 k \ T A (k, M B ) T A (k, M B ) 



dP^ J (2-nfE^k) \ ml - (P - kf) 2 J 

fk r A (k, M B ) 2(P - kfr T A (k, M B ) 

, 2 



This is a familiar result, which will be generalized to the spin 1/2 case later. 



(3.10) 



E. Symmetrized two channel equation for equal mass scalar quarks 

If the quarks have equal mass (mi = m 2 = m), and the bound state mass is positive 
and not too small, a symmetrized two channel equation is needed. The two channels will 
be labeled 1 and 2 depending on whether the quark or antiquark is on mass-shell, and the 
symbol pf denotes that the particle is on its positive energy mass-shell, (ie. pf 2 = m 2 and 
Pq = E(p) = \Jm 2 + p 2 ). Starting from Eq. ( |3.7| ), and suppressing the subscript A, the 
vertex functions for the two channels are denoted 
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ri( P ,M B ) = r(p+ P2 ), r 2 (p,M B ) = r( Pl ,pt 



(3.11) 



With this notation, the symmetrized two channel equation for equal mass scalar "quarks" 
with a confining interaction can be written 



Up, M B ) 



-m 



d 3 kV i:j (p,k) 



, (2tt) 3 Ej(k) 
2mCT i (p,M B ) 



[P - k+) 2 



P 



Pi 



m 2 — (P — p- 

where % and j label which of the two quarks is on-shell, and 

k+ = {E(k),(-y^k} 



(3.12) 



(3.13) 



is the momentum of the on-shell quark. Note that the strength of the V^- term has been 
multiplied by 1/2, reflecting the fact that the interaction is an average of the strengths 
in two channels which are equal in the nonrelativistic limit. This equation uses the same 
subtraction for both the i = j and the % ^ j terms. This prescription differs from that 
previously used in Ref. [1[]. In this work the kernel below will not, in general, be singular 
when i ^ j, and the subtraction used above is sufficient to preserve the nonrelativistic limit 
(see below). 

In order to complete the description we need to specify the form of covariant interaction 
Vij. A natural choice that reduces to the correct nonrelativistic limit is [|l| 



Vij(p, k) = V A {qi 



8na 



(3.14) 



where the four-momentum transfer depends on whether or not i = j: 



q 2 n = q 2 22 = (E(k)-E(p)) 2 -(k-p) 
q 2 12 = q 2 21 = (M B -E(k)-E(p)) 2 -( 



A similar form could be used for the kernel V$ (which we will not need) 



Vijip, k) = V s (qi 



— 57TO- 



+ 



4e 2 



- ^ 2 ) 2 (qI - £ 



2\3 



(3.15) 



(3.16) 



However, the form ( ^.14[ ) has two drawbacks. First, at large p ~ k the kernel converges 
slowly, and the equation is ultraviolet divergent. In Ref. U a form factor was introduced to 
regularize this divergence. Second, using this form it is difficult to regularize the infrared 
(q 2 = 0) singularities that appear in the e = limit. In the nonrelativistic case the infrared 
singularity occurs only at q = and can be regulated by the 5 function subtraction in 
Eq. (|2.3| ). However, in the relativistic case infrared singularities occur not only when q^ = 0, 
but also (for the % ^ j kernels) when the momentum transfer is light-like, so that q 2 = but 
g M 7^ 0. These " off- diagonal" singularities are not regulated by the subtraction term, and 
their removal spoils the simplicity of this approach IJ. 
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Since the role of Va is to model the linear interaction, and the principle requirement is 
that it reduce to the correct nonrelativistic limit, both of these problems are eliminated very 
simply if Va is defined as follows 

VAiXli) = -«j + (p-«,)vj» (3 ' 17) 

where P is the total four-momentum of the bound state. This form has the following 

advantages: 

(i) the denominator is not singular unless both q 2 and P ■ q are zero, so the singularities 
are restricted to q^ = 0; 

(ii) no ultraviolet regularization is needed; 

(iii) the interaction does not depend on the bound state momentum P in the bound state 
rest frame; and 

(iv) it has the correct nonrelativistic dependence on q 2 . 



One disadvantage of the form ( [3.171) is its dependence on the total momentum P of the 
particle pair. However, since since this kernel confines particles in pairs that can not be 
separated, they are naturally associated as a pair and we do not view this as a serious 
limitation. Another feature of the form ( |3.17| ) is that its off-diagonal couplings are singular 
only when W = 2E{p) (because k + p = also). This is only possible for excited states and, 
as we will prove below, confinement requires the vertex function to be zero at this point, 
controlling this singularity automatically. 

The introduction of the definition ( |3.17|) considerably simplifies the solution of the rela- 
tivistic equations ( 3.12Q , but will introduce electromagnetic interaction currents if the photon 



four-momentum is not zero. These will be discussed in a subsequent paper. 

Both Eqs. ( |3.7| ) and ( |3.12| ) have the correct nonrelativistic limit with confinement. Con- 
sider the one-channel Eq. (|3.7|) first, and let rri\ and m 2 — > oo. Then the energy trans- 
ferred by the on-shell quark, E\(k) — E\(jp) — ► and Va{q\.i) — > V^q). Furthermore, if 
Mb = vn,2 + mi + E, then to first order in the small quantities k 2 and turE, the relativistic 
propagator reduces to 

1 mR (3.18) 



m\ — k\ m 2 (k 2 — 2m rE) ' 

and substituting this into Eq. (|3.7| ) gives Eq. fl2.26| ). In the two channel case qu — > q\% 
as m — > oo and the kernels Vn — > Vu- Since the subtraction in the two channels is also 
identical, the contributions from the two channels are equal and the coupled equations reduce 
to the single Eq. (|2.26|) . 



F. Proof of confinement 

While one can visualize the potential in the nonrelativistic case and get a picture of the 
physics, it is less possible to visualize the covariant interaction. What are the criteria with 
which one can judge whether a given interaction really confines'? If the particles are bound in 
a state of total mass larger than the sum of the masses of the constituents (M B > m 1 + m 2 ), 
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FIG. 5. The confinement condition for the Gross vertex function. 



the bound state could in principle decay into free constituents. Confinement prevents this 
from happening in one of two possible ways: (i) the quark propagators will not have any 
physical mass poles |H| , or, as we will now prove for this model, (ii) the vertex function will 
vanish when the quarks are simultaneously on-shell. 

The proof is identical to the nonrelativistic proof given above and we will summarize it 
only for the one channel equation. Setting C — 0, the one channel bound state Eq. 
can be written 



T A (p, M B ) = -2m 1 m 2 J 



d 3 k V A (p, k) 



r A (k,M B )-T A (p,M_ 



B, 



in- 



/to 



+2m 1 m 2 T A (p, M B ] 



d 3 k V A (p, k) 



p\-k\ 



(m: 



(ml 



(3.19) 



Since the first quark is on-shell, the second quark is on its positive energy mass shell when 
the magnitude of the relative three-momentum |p| = po is 



M, 



B 



(3.20) 



This occurs when is given by 



am 2 b pI 



Ml 



(mi + m 2 ) 



Ml 



(mi - m 2 f 



(3.21) 



As in the nonrelativistic case, the singularity at p = k is integrable, and hence the second 
term on r.h.s. of Eq. ( |3.19 ) will be singular at p = p p = p (where p is a unit vector in 
the direction of p) unless 



r( Po ,M B ) = 



(3.22) 



Therefore, the vertex function vanishes when both particles are on their mass shell. This 
condition is illustrated diagrammatically in Fig. || 

Note that the subtraction term in Eqs. ( |3.12| ) and ( |3.19| ) plays two central roles: (i) it 
regularizes the singular interaction at p = k and and makes it zero at r = 0, and (ii) it is 
singular when p\ —>■ r??|, forcing condition ( |3.22 ). The subtraction term is essential to the 
self consistent description of confinement. As in the nonrelativistic case the proof did not 
depend on the specific form of the interaction. 

We now discuss how confinement affects the stability of bound states under external 
disturbances. 
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FIG. 6. The Born term, which can not exist if the quarks are confined. 




G. Excitation of bound states 



A consistent description of confinement implies that two free quarks can not be liberated 
from a bound state, even under the influence of an energetic external photon or other probe. 
This requirement implies that the usual Born term (shown in Fig. |6]) is either cancelled by 
the rescattering term, or is a diagram that does not exist in the formalism. If the Born term 
does not exist, the rescattering term, illustrated in Fig |7|, must be zero if the final state 
quarks are all on-shell. How are these restrictions built into the formalism? 

When particles are confined there are no free two-particle states and the two-body propa- 
gator must always include an infinite number of interactions. Since there are no free particle 
states, a perturbation theory for confined particles built around the free propagator can not 
be constructed. This feature is built-in automatically if the two body propagators satisfy 
homogeneous integral equations with no free particle contribution. 

To illustrate these ideas we review the formalism for the scattering amplitude, and its 
relation to the two-body propagator. It is convenient to work with the scattering form of 
the equation. In operator notation, Eq. ( |2.17| ) is: 



where Go(k, k', P) is the free two body propagator [containing a factor of S 3 (k — k')], inte- 
gration over d 3 k and d 3 k' is implied, and we have dropped the subscript S for simplicity. 
The parameter rj was introduced in the discussion following Eq. fl2.19|) and is very small, 
approaching zero as e — > 0. 

Now the dressed propagator, G is related to the scattering amplitude M by 



M(p, p', P) = r, V(p, p') - V(p, fc)G (k, k', 



P)M(k',p',P) 




(3.23) 



G(p, p', P) = C G (p, P, P) - G (p, k, P)M(k, k', P)G (k', p', P) , 



(3.24) 
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where £, to be determined, is a parameter proportional to the strength of the free particle 
scattering. If the potential confines there should be no inhomogeneous term and ( = 0. To 
determine ( and the equation for G, substitute ( |3.23| ) into ( |3.24j ) giving 

G(p, p', P) = C G (p, P, P) - ^ (p, k, P)V{k, fc')G (k', p', P) 

+G ( P , k, P)V(A;, fc')G (k', k", P)M(k", k w , P)G (k w , p', P) 
= C G (p, p', P) + (C - 7/) G (p, k, P) fc')G (k', p', P) 

-G (p, k, P)V(A;, fc')G(k', p', P) . (3.25) 

The second term is eliminated by choosing ( = rj, and gives familiar equations for the dressed 
propagator 

G(p, p\ P)= V G (p, P ', P) - G (p, k, P)V(k, k')G(k', p', P) 

= 7] G (p, p', P) - G(p, k, P)y(A;, fc')G (k', P ', P) , (3.26) 

where the second form parallels the second form of Eq. ( |3.23| ). 

The interpretation of equations (|3.24 ) and (|3.26j ) for the dressed propagator follows from 
the interpretation of Eq. ( |3.23| ) for the scattering amplitude. As e — > 0, the parameter rj — >• 
and the inhomogeneous term vanishes. In this limit both the scattering amplitude and the 
propagator satisfy homogeneous equations. 

The inelastic scattering amplitude can be obtained from the dressed propagator by strip- 
ing off the final free propagators, and is 

J(p, P, q) = Go '(p, k, P + g)G(k, p', P + q)J(P + q, P)^(P) 

= {rj + M(p, k, P + g)G (k, p', P + q)} J(P + q, P)*(P) . (3.27) 

Here the first term proportional to r\ is the Born term shown in Fig. |||, and we see that there 
is no Born term in the limit of exact confinement (ie. 77 = 0). Furthermore, in the presence 
of confinement the scattering matrix satisfies the same homogeneous equation satisfied by 
the bound states [Eq. ( p.23| ) with 77 = 0], and an extension of the proof given in Subsec. F 
above shows that the scattering matrix in Fig. [J must be zero if both final state quarks are 
on shell. 

We have constructed a self-consistent description of confinement within the context of 
relativistic field theory. 



H. Generalization to Fermions 

If the quarks have spin, the kernel in the spectator equation will be an operator in the 
Dirac space of the two quarks. This operator can be written 

3 

V(p, k) = £ on O a O l2 Viip, k) , (3.28) 

i=l 

where the Dirac matrices O, which operate on the Dirac indices of particles 1 and 2, describe 
the spin dependent structure of quark- ant iquark interaction. The oti are parameters deter- 
mined either empirically (by fitting the spectrum), from lattice calculations, or from the 
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i=l 

a. 



theory. In this paper we consider only three possible spin structures: scalar 0\j = ILj, pseu- 
doscalar 2 j — J5j, an d vector O^j = 7/y/2. With this notation the one channel spectator 
equation for spin 1/2 particles with constant masses mi » m 2 is given by 

r(p ' p) = - / wfm ? Vi{p ' k)0a(mi + M {^ri} ^ - *>°» (3 - 29) 

where the quark has mass m 1 and is on shell, so that kf 2 = m\ = pf 2 , and the antiquark 
has mass m 2 . Therefore, the momentum transfered by the interaction is 

(p+ - fc+) 2 = (^(p) - ^(fc)) 2 - (p - k) 2 = q. (3.30) 

As in the nonrelativistic case, we consider a kernel composed of linear, constant, and one 
gluon exchange (OGE) pieces. The interaction kernel for the linear part of the potential, 
V z , is 

3 

Vz(p, k) =Y^ ® Li O a O i2 V L (p, k), 

3 lLilL 2 + ^751752 + - a v j^ij^) V L (p, k) , (3.31) 
where Vl(j>, k) is 

V L (p, k) = V A (q n (p, k)) - Ex(k) 5 3 (p -k)J d'k'^^pi . (3.32) 

In this work we employ a pure scalar linear interaction, a s = 1, a ps = a v = 0, but in later 
calculations the coefficients aj will be determined empirically. The one gluon exchange and 
constant interactions will be pure vector 

Vc(q)=l»il%C } (3.33) 

where 

vr(ff) = -or - ^) %) 

= j 16?r2 / 3 (3 34) 

l ^ g 2 V-A 2 ln(r + | g 2 |/A 2 QOT )' ^ } 

where d = 12/(33 - 2JV» = 12/27, the color factor of 4/3 has been included, A = 1 GeV, r 
= 2, and Xqcd = 200 MeV. In previous work [0 quark propagators with constant masses 
were used. In this work we parametrize the quark propagator by 

where m{p) is a mass function for the quark, to be defined later. 
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FIG. 8. Propagator poles in the complex ko plane. 



If the constituents are identical or close in mass and the equations are to be applied to 
the description of nearly massless bound states, the four channel equation should be used. 
Numerical solutions the four channel equation will be presented in this work. 

The four channels are defined by the constraints in the four-momenta k\ and k 2 arising 
from the requirement that both the quark and the antiquark be constrained to both their 
positive and negative energy mass-shells. A formal way to obtain the equations is to integrate 
over the internal energy ko by averaging the contributions from the quark and antiquark poles 
in both the upper and lower half ko complex plane, as illustrated in Fig. |[ This averaging 
is needed to ensure charge conjugation (particle-antiparticle) symmetry, and leads to four 
coupled equations. However, even though the form of the equations is obtained in this 
way, we emphasize that the equations are theoretically justified by the argument that the 
singularities in the interaction kernel omitted in this procedure tend to be cancelled by 
other higher order terms which would otherwise have been neglected, and that this leads 
to covariant equations with the correct nonrelativistic limit. The inclusion of the negative 
energy poles, neglected in other applications of the symmetrized equations ||, is required 
in cases where P — ► M. 

The four constraints are conveniently identified by the notation 

k] = {sE(k),(-) j+1 k}. (3.36) 

which generalizes that introduced in Eq. ( |3.13| ). Here the superscript s = ± denotes either 
the positive or negative energy mass shell constraints. Then, introducing the projection 
operators 

A(k) = m(k) + ^, (3.37) 
and defining the four channel vertex functions 

r 1 (p,M B ) = T(p s 1 ,p 2 ) 

T s 2 (p,M B ) = T( Pl ,p s 2 ), (3.38) 



and wave functions 



A(pf)rf(p,M B )A(-p2) 
m(p 2 ) 2 ~ v\ 



*-r« m \ A( Pl )r-( P ,M B )A(-^) 

* 2(p ' Mb) = m( Pl y- P j ' (3 - 39) 
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permits us to write the four-channel spectator equation in the following compact form 



d 3 k 



2- J (2^2E(k){ V ^- k nn^M B )-m P ,M B ) 



J I' 



- 25 ij 5 sr Vf{ P - k^Mk, M B ) lv - C Tm ^(p, M b )Y 



(3.40) 



where the r.h.s. of the equation now sums over both positive and negative energy con- 
tributions (r = ±) from each quark (j = ±). The Kronecker Sij5 sr functions restrict the 
one gluon exchange interaction to the diagonal channels (where the same particle is on the 
same mass shell before and after the interaction). Inclusion of the one gluon exchange in 
off-diagonal channels leads to numerical instabilities, which in principle can be handled by 
using more grid points in numerical integrations. Restricting this interaction to diagonal 
channels eliminates these singularities from the gluon propagator. 



I. Charge conjugation invariance 



The final task is to show that Eq. fl3.40|) is invariant under the charge conjugation oper- 
ation 

T c ( Pl ,p 2 ) =CT T (p 2 , Pl )C- 1 . (3.41) 

This is done by proving that both T and T c satisfy the same equation. 

First note that, when particle 1 is on shell, interchange of p\ and p 2 gives 

T{(p,M B ) =T(plp2) ^r{p 2 ,p{) = r s 2 {-p,M B ) (3.42) 

and is equivalent to 1 «-> 2 and p — > — p. Then 

^ c (p, M B ) = C% T (-p, M B )C- 1 

^ c (p, M B ) = C^ T (-P, M B )C- 1 . (3.43) 

Finally, the Dirac direct products 11 <g) 11, 7 M ® 7 M , and 75 <g) 75 are invariant under C. 
Hence, changing k — > — k and performing the transformations ( |3.41| ) and ( |3.42| ), shows that 
Eq. ( p. 40 ) is also invariant. Therefore the charge conjugation eigenstates, labeled by rj = ± 

r*( P , M B ) = r;(p, m b ) + v rf ( P , m b ) , (3.44) 

are solutions of the equation and charge conjugation symmetry is proved. 



J. Dynamical quark mass 

The dynamical quark mass function is the solution of the Dyson-Schwinger equation. In 
NJL-type models, this one-body equation for the spontaneous generation of quark mass and 
the two-body bound state equation for a state of zero mass become identical in the chiral limit 
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(when the bare quark mass is zero). In this limit the quark mass function and the bound 
state wavefunction for a massless pseudoscalar bound state are identical, and spontaneous 
symmetry breaking assures the existance of a massless pseudoscalar bound state. 

In this paper we adopt a slightly different approach. We will first choose a convenient 
mass function, and then require that the two-body equation for a massles pseudoscalar bound 
state automatically have a solution when the bare quark mass is zero. In this case the quark 
mass function and the wave function for the massless Goldstone boson will not be identical, 
but at least the existence of the Goldstone boson in the chiral limit is assured. We will 
define the quark mass function of flavor / by 

mf(p) = m jr + c(m°j-)f(p), (3.45) 

where is the current quark mass of flavor /, and f(p) is a universal function defined by 

/(p) s HTA^' (3 ' 46) 

The function c(m°) can be thought of as a polynomial in powers of m". This is the typical 
structure of the mass function which is usually obtained from the solution of the one body 
equation. 

The reason for not solving the one body equation, in our case, is two fold. The first 
problem is the difficulty of incorporating one gluon exchange into the one body equation. 
Because of the on-shell constraint in the loop momenta, the one gluon exchange interaction 
leads to an ultraviolet divergence. The second problem is associated with our choice of 
infrared regularization of the linear interaction. The infrared singularities are regulated by 
the P ■ q term in the denominator of the linear interaction Eq. (13.17 ), and this would imply 



that the resultant mass function is a function of two arguments, i.e. m = m(p 2 , p 2 ). This is 
unacceptable, and rather than forsaking important features of the model such as confinement 
and asymptotic freedom, we choose to model the quark mass functions. 

The form ( 3.45|) guarantees that at large momenta, quark masses go to their current 



quark mass values as dictated by asymptotic freedom. In the chiral limit the quark mass 
reduces to 

m x (p)=c(0)f(p) (3.47) 

We fix the constant c(0) by requiring that the pion bound state equation, using the mass 
function Q3.47 ), give a massles solution. This insures that a massless pion exists in the chiral 



limit when mj = 0. Next we choose a value for the light current quark mass, m° u = mj, and 
fix c(m°) so that the two-body equation gives the correct value for the physical pion mass. 
This also fixes the value of the on-shell quark mass away from the chiral limit. Similarly, we 
choose mP s and fix c(m°) by fitting the kaon mass. For three flavors it is therefore sufficient to 
have a function c{m°f) which is a polynomial of order 2 in m". As new flavors are introduced 
the order of the polynomial accordingly can be increased. 

To summarize, we have 6 mass parameters: m°, m", c(0), c(m°), c(m° s ), and A. In practice 
we fix A at one GeV and choose the current quark masses m° ul and m° s to be near the values 
expected by current theory. We then adjust the c's to give the a zero mass pion in the 
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chiral limit, and a real pion and kaon with the observed masses. This process is repeated 
for different values of the current quark masses and the potential parameters a and C until 
satisfactory values for the constituent quark masses and the spectrum of excited pions is 
obtained. The final values of the parameters will be given in the next section. 

Having outlined the features of the model, we now turn our attention to the details of 
the pseudoscalar bound state equation with spin. 



IV. PSEUDOSCALAR CHANNEL 

The bound state vertex function has the following structure 

X = Xcolor ® Xflavor ® Xspin- (4.1) 

The color space vertex function is a Kronecker delta function, 5 c d, which reflects the color 
singlet nature of the bound state. The flavor space vertex function is the matrix in 577(3) 
matrix space, which chooses the right flavor combination of the meson under consideration. 
Indices /, g refer to up down and strange quark entries (it, d, s — 1, 2, 3) of A\ For example, 
[A + ]«d = [^ + ]i2- For a general meson type i, the bound state vertex function is 

xU,fg,cd( k u ki) = \* fg r a/3 (k u k 2 ). (4.2) 

where a, (3 are Dirac indices (to be suppressed in the following discussion). The most general 
form for the spin-space part of the vertex function for pseudoscalar mesons is 

r(h, k 2 ) = 75 {r + f r x + 1 r 2 + [ % f] r 3 } , (4.3) 

where Tj = rj(fei,A;2) are scalar functions. The dominant contribution to the bound state 
vertex function comes from the first term of (14.31), 



T{k u k 2 )~l^»{k u k 2 ), (4-4) 

This approximation, which is exact in the chiral limit when P = and ni\ = m 2 , will be 
used for the pion and kaon bound states in this work. 

Assuming ( |4.4| ), multiplying the four channel equations for pseudoscalar mesons by 75, 
and taking the trace, gives the following approximate coupled equations for pseudoscalar 
states 

1 r d z k ( r 
r?(p, M B ) = - - E J (2n)32E . (k) \ V iI(Pi k ) [mW^ M b) - FMTXp, M b 



+ 65 tJ 6 sr V g (p - A;)F i (^)r j r (k, M B ) + 2F i (pJ)CT?(p, M B ) , (4.5) 



where the four channel wave functions r|(p, Mb) are obtained from r as shown in Eq. ( |3.38| ) 
and 
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mi) 
mi) 



m 1 m 2 (h) + k[- k 2 

ml{k 2 ) - kl 
mi(fei)m 2 + k\ • k T % 



m\{k 2 ) - k\ 



(4.6) 



where mj(fc[) = m.j(— hi-) = mi. For future reference we record the four-momentum q^J 
(pi — ki)lj exchanged between the two quarks. This depends on the initial and final channel. 
The distinct cases are: 



q[ s 1 =q 22 r '- S = (rE(p)-sE(k), p - k) 
q[°= {rE{p) + sE(k) — Mb, p — k) 
q%= (M B -rE(p)-sE(k), p - k) 



(4.7) 



The solution of Eqs. (|4.5|) for a realistic choice of the parameters will be discussed in the 
next section. 

Before turning to this discussion, look at the coupled equations in the chiral limit, when 
P = and the dynamical quark masses are equal, so that m\{k) = m 2 (k) = m(k). In this 
limit, k\ = —k 2 , and expanding to order P ■ k\ gives 



_ m{k[)m{P-kl)+kl-P-kf _ 1 - 2mm' _ 1 _ 

rn?{P-kl)-{P-klY 2 -4mm' 2 21 2j 



(4.8) 



where m = m(±fc?) and m' = dm(±k) / dk 2 \^ 2 =m?)- Hence, using charge conjugation sym- 
metry fl3.4lQ , the four coupled equations Q4.5|) reduce to only two equations in the chiral 
limit. These coupled equations are 



r+( P ,o)= - 



d 3 k 



v+(p,k) [r+(k,o)-r+( P ,o) 



r;( P ,o) 



(2tt) 3 2E(k) 

+ v^p, k) [r-(k, o) - r+( P , o)] + qv 9 ( p - fc)r+(k, o)| + 2cr+( P , o) 

r d 3 k { 
-/(^2l(I){^- fc) ^( k . )- r x-(P.O)] 

+ V4 P , k) [rj(k, o) - r-( P , o)l + 6V S ( P - k)r-(k, o) 1 + 2cr-( P , o) , (4.9) 



where 



V±(jp,k) 



8trj 



(p-k)* + (E(p)=F 
Note that these two equations are symmetric under the interchange 

and hence reduce to one equation for r x = Ft. = ±r 



(4.10) 



(4.11) 



x 



x 
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FIG. 9. Quark mass functions mjip) = M(p 2 ) are shown for up/down, and strange quarks. 
On-shell quark masses are m u ^ = 360 MeV, and m s = 588 MeV. At large momenta quark mass 
values approach to d = 5 MeV, and ra" = 100 MeV. 



r x ( P ,o)= - J 



d 3 k 



[V + (p, k) ± 7_(p, k) + 6V g (p - k)) T x (k, 0) 



2nf2E(k) 

[V + (p, k) + V-(p, k)] T x (p, 0)) + 2CT x (p, 0) 



(4.12) 



where the sign of the V_ term depends on the sign in the relation (|4.11|) . Since the tt is 
even under charge conjugation symmetry, the plus sign is the correct one to use. 
Recalling Eq. (|3.47| ), the energies E in Eq. (|4.12|) depend on the constant c(0) 



E(p) = Jc(oyp(p) + P * 



(4.13) 



and this is adjusted to insure that the Eq. (|4.12|) has a solution. Once c(0) has been fixed, 
Eqs. ( |4.5|) are solved for various values of the bare quark masses and the "mass functions" 
c(mj), and all parameters are adjusted to give a reasonable spectrum. 

Having outlined the features of the model we next present the results for mass functions 
of quarks and vertex functions for bound states. 



V. RESULTS 

The quark mass functions are shown in Fig. |9|. The on-shell quark masses rrif are given 
in Table I. At large momenta, the quark mass values approach the bare quark masses m' 
shown in Table II. The other mass parameters and bound state parameters are also shown 
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TABLE I. Summary of results 



Observable Calculated Experimental 

140 MeV 139.6 MeV 

m v 320 MeV 

m* 1118 MeV 1300± 100 MeV 

m K 495 MeV 495 MeV 

m x 376 MeV 

m u = m,d 360 MeV 

m, 588 MeV 



in Table II. The parameter A which determines the scale of mass function was fixed at A = 1 
GeV and not adjusted during the fits. The third line in Fig. ||| is the momentum p, and the 
intersection of this line with the quark mass function gives the constituent quark mass. 



In Figs. [10] and [TT] the ground and first excited state vertex functions of the pion are 
shown. Here we show the vertex functions as a function of the variable p s j = sE(p) = spo. 
Note that p is positive for positive energy states (s = +) and negative for negative energy 
states (s = — ). Because of the symmetrization, the positive energy quark vertex function 
is the same as the negative energy anti-quark vertex function up to an overall phase (+ 
for states even under charge conjugation and — for odd states). Also note that the curves 
are not continuous because the argument p$ can not take values between (— m, +m). In 
Fig. [TJ we present the excited state vertex functions on a logarithmic scale. The location of 
the first node is exactly where both quarks are simultaneously on shell. Therefore, although 
kinematically allowed, the excited state of the pion can not decay into a free quark- ant iquark 
pair. This numerical result is a consequence of the confinement condition ( j?. 23J . 

In Fig. [13| we present the non strange-eta (the isospin zero uu + dd combination) ground 
state vertex functions. Note that these are odd under charge conjugation. The kaon vertex 
functions are shown in Fig. |14[ Since the kaon is formed from a quark and antiquark of 
unequal masses, the particle-antiparticle symmetry is lost and the negative and positive 
energy solutions have a differerent shape and size. 

The mass function and the pion wave function in the chiral limit are shown in Figs. |15 
and |I6. 



VI. CONCLUSION 

We have shown that a relativistic generalization of the Schrodinger equation with linear 
interaction leads to the Gross equation. It is not possible to write a Bethe Salpeter equation 
that gives the correct linear interaction in the nonrelativistic limit. We have proved that the 
relativistic generalization of the linear interaction leads to vanishing vertex amplitudes when 
both of the constituents are on-shell. This guarantees that the bound state does not decay 
to its constituents. This mechanism of confinement follows from insisting on the correct 
nonrelativistic limit. The model incorporates asymptotic freedom through the inclusion of a 
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TABLE II. Values of the parameters 



Parameter Value 



m° 100 MeV 



ml 5 MeV 



s 

c(0) 0.429 (Gev) 3 

c(m°) 0.400 (GeV) 3 

c(m°) 0.657 (GeV) 3 
a 0.4 (GeV) 2 

C 0.4929 
A 1 GeV 




FIG. 10. The four-channel vertex functions for the ground state of the pion. 
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FIG. 11. The four-channel vertex functions for the first excited state of the pion. 
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FIG. 12. The two positive energy vertex functions for the first excited state of the pion. The 
second node is due to the excited state, and the first node assures that the bound state does not 
decay. 
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FIG. 13. The four-channel vertex functions for the non-strange ij. 
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FIG. 14. The four-channel vertex functions for the ground state of the kaon. 



27 



0.5 




1e-03 1e-02 1e-01 1e+00 1e+01 1e+02 

p 2 (GeV 2 ) 

FIG. 15. The chiral limit of the quark mass function M(p 2 ) = m x (p). The on-shell quark mass 
is m x = 376 MeV. At large momenta quark mass function approaches 0. 
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FIG. 16. The chiral limit of the pion ground state vertex function. 
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vector one gluon exchange interaction, and quark mass functions that approach the current 
quark values at infinite momentum. There are no cut-off's or ad-hoc form factors involved, 
and the linear interaction involves only one coupling parameter. The approach give a good 
description of the pion, kaon, and eta. 

It remains to use this formalism to describe the full meson spectrum. 
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APPENDIX A: NUMERICAL METHODS 



Solutions of integral equations are performed by first discretizing the integrals 

„ n 



i=i 



where Wi are integration weights for grid points gj. In order to map the grid points and 
weights from interval (—1, 1) to (0, oo) we use the arctangent mapping (Ref. [pT| .|l2|1 ) 



y(x) = R r , 



71 



R d tan ( -(1 + x] 



1 + 



R,i 



Rmax Rmin 



tan + x )^j 



(A2) 



where 



It follows that 



R,i 



Rmed Rmin 



R 



max Rmed 



(R, 



max ± ^mtn 



R„ 



y(-i) 



Rn 



2/(0) = R 



med 



y{i) = R, 



(A3) 



(A4) 



Therefore, one can safely control the range (R m in,Rmax) and distribution (R m ed) of grid 
points. With this discretization procedure, continuous integral equations are transformed 
into nonsingular matrix equations. 

The spectator equation is an eigenvalue problem, where the eigenvalue is the mass of the 
bound state. The equation can be brought into the following form 



N 



'£[HM(pi,P j )-l]*(p j ) = 0. 



(A5) 



where M is the bound state mass, and Pi (i = 1...N) are grid points in momentum space. 
Therefore, H is an N x N matrix and $ is a vector of dimension N, which leads to the 
following matrix equation 
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[H M - 1]$ = 0, (A6) 

where M is unknown. Start by making an initial guess for M. In order to find the ground 
state, one should start with an initial guess near the expected value of the ground state 
mass. The next step is to see whether the initial guess leads to a consistent solution. The 
most efficient way of checking whether a given matrix has a specific eigenvalue is through 



the method of inverse iteration, as suggested in Refs. |TT],|T2[. First construct an arbitrary 
vector x° 



where $j, i = 1..N, satisfy 



N 

X° = E C ^ (A7) 



[H M -u i ]$ i = Q, (A8) 



where uii, i = 1..N are eigenvalues of the Hm matrix. It should be emphasized that eigen- 
values which are not equal to 1 have no physical meaning, for they do not correspond to a 
solution of the equation (Eq. |A~6| ). Next, construct 

K = 7T-^- (A9) 

Letting K operate on state x° n times produces 

When the number of iterations n is sufficiently large (usually around ten), the dominant 
contribution to % n comes from the eigenvector whose eigenvalue ujj satisfies \ujj — 1| < 
\uji — 1| for all i — 1 • • -j — 1, j + 1 • • ■ N. Therefore, 

X n+1 ^^—rX n - (All) 
Using the eigenvector x™, which is proportional to $j, the eigenvalue Uj can be found from 

n\ tt n 

w , = x + x . (A12) 

X X n 

If ujj is close enough to 1, then one has a self consistent solution. This method has the benefit 
of directly singling out the eigenvalue closest to the initial guess, rather than finding the 
largest eigenvalue as in the case of straight forward iteration. Excited states can similarly 
be found by varying the initial guess M towards higher values. There is only one matrix 
inversion involved. Distribution of the grid points in momentum space is done by the 
arctangent mapping. The typical number of momentum space grid points used in order to 
obtain stable solutions is around 40. 
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